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Abstract 

Thermal conductivity of superconducting MgB2 was studied in 
both the superconducting and the normal state region. The latter 
is almost equally determined by the electronic - and the lattice con- 
tribution to the total thermal conductivity. In the superconducting 
state, however, the lattice contribution is larger. The electronic ther- 
mal conductivity below was derived from the experimental data 
considering the Bardeen-Rickayzen-Tewordt theory together with the 
model of Geilikman. The analysis shows that electron scattering on 
static imperfections dominates. 

key words: MgB2 superconductivity, thermal conductivity, spe- 
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1 Introduction 



Superconductivity with a remarkably high transition temperature Tc ~ 39 K 
was recently discovered in MgB2 The subsequent investigation of the 
boron isotope effect by Bud'ko et al. p[ revealed a partial isotope exponent 
q;b ~ 0.26 (corresponding to ATc = 1 K) which appears to be consistent with 
a phonon-mediated BCS superconducting (SC) mechanism. Also other ex- 
perimental reports as e.g. on specific heat [§, Q argued that their data can be 
accounted for by a conventional, s-wave type BCS-model. Investigations of 
the SC gap of MgB2 by means of Raman scattering are also consistent with an 
isotropic s-wave gap with a moderate coupling 2 A ^ A.lksTc 0. However, a 
theoretical analysis of the temperature dependence of the upper critical field 
Hc2(T) in terms of Eliashberg type models by Shulga et al. demonstrates 
that the shape and the magnitude of the upper critical field of MgB2 can 
definitely not be accounted for by an isotropic single band model, but may 
successfully be described within a multi-band Eliashberg model with various 
options. A careful calorimetric investigation of the SC parameters of MgB2 
by Wang et al. [j^ gave even more direct evidence against the arguments 
for simple isotropic BCS type superconductivity: "The nearly quadratic de- 
pendence of C(T) versus T at T ^ Tc, its non-linear field dependence, and 
the discrepancy between the electron-phonon coupling constant X^p as deter- 
mined by the renormalization of the electron density of states (Aep ~ 0.6) 
and by McMillan's equation for isotropic superconductors (Aep ~ 1.1), are 
inconsistent with a single isotropic gap" . Direct hints for a non BCS temper- 
ature dependence of the gap energy A(T) were also obtained by tunneling 
experiments on MgB2/Ag and MgB2/In junctions 

Thermal conductivity A is one of those transport coefficients which ex- 
hibits non- zero values in both the normal and the SC state. The temperature 
dependence of A allows to distinguish between the most important interac- 
tions present in a superconductor. In particular, the interaction of electrons 
with phonons are recorded in the magnitude of A(T). Moreover, scatter- 
ing of these particles by static imperfections like impurities, defects or grain 
boundaries are reflected. 

The aim of the present work is to derive the temperature dependent 
thermal conductivity of MgB2 and to analyse the data with respect to the 
electronic - and the lattice thermal conductivity both in the normal and the 
SC state. Moreover, we present resistivity and specific heat measurements 
in order to characterise the investigated sample. 
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2 Experimental 



A MgB2 sample of about 1.3 g was synthesized by direct reaction of the el- 
ements. The starting materials were elemental magnesium (rod 99.9 mass 
% nominal purity) and boron (99.5 % powder, crystalline, < 60 mesh, 99.5 
mass %). The elements in a stoichiometric ratio were enclosed in a cylindric 
tantalum crucible sealed by arc welding under argon atmosphere. The tan- 
talum crucible was then sealed in an iron cylinder and heated for one hour 
at 800° C and two hours at 950° C in a furnace. The sample characterized 
by x-ray diffraction show pure MgB2 phase; only one very weak peak due to 
extra phase has been found. 

The thermal conductivity measurement was performed in a flow cryo- 
stat on a cuboid-shaped sample (length: about 1 cm, cross-section: about 
2.5 mm^), which was kept cold by anchoring one end of the sample onto 
a thick copper panel mounted on the heat exchanger of the cryostat. The 
temperature difference along the sample, established by electrical heating, 
was determined by means of a differential thermocouple (Au -|- 0.07 % 
Fe/Chromel). The measurement was performed under high vacuum and 3 
shields mounted around the sample reduced the heat losses due to radiation 
at finite temperatures. The innermost of these shields is kept on the temper- 
ature of the sample via an extra heater maintained by a second temperature 
controller. 

Resistivity data were taken from bar shaped samples applying a standard 
4-probe d.c. technique at temperatures down to 0.5 K and in magnetic fields 
up to 12 T. 

Specific heat measurements were carried out on a sample of about 1 g in 
the temperature range 5 K-50 K using a quasi adiabatic step heating tech- 
nique. 

3 Results and discussion 

In order to give a direct proof of the SC bulk properties of our MgB2 sample 
prior to the transport measurements we checked the specific heat . These 
specific heat measurements performed in zero-field, 1 and 9 T showed reason- 
able agreement with results previously reported 0, 0, 0] . The thermodynamic 
mean transition temperature Tc of our sample is 37.5 K. As already noted in 
the introduction there have been distinctly different and partly controversial 
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conclusions suggested in Refs. ^, 0, although their raw data are in fair 
agreement with each other. Therefore, we show in figure |l| the temperature 
dependence of the electronic specific heat, Cei(T) versus T, obtained by sub- 
tracting the lattice heat capacity deduced from the 9 T specific heat data. Of 
course, 9 T are insufficient to obtain a complete suppression of superconduc- 
tivity in MgB2, but there is already a large reduction of the order parameter 
combined with a dramatic broadening of the transition (see Ref. 0] for com- 
parison of 10, 14 and 15 T data). Thus, we obtained a Sommerfeld coefficient 
of the normal-state electronic specific heat 7 ~ 2.4(2) mJ/molK^ by extrap- 
olating the 9T data in C/T versus from 30 - 100 to zero temperature 
which is already close to 7 = 2. 7 ±0.15 mJ/molK^ obtained from the 14 and 
16 T measurement by Wang et al. The important point to emphasize in 
figure |l| is the non BCS-like temperature dependence of the SC-state elec- 
tronic specific heat. In fact, we observed a similar deviation from a simple 
BCS temperature dependence (solid line, figure |l]) as previously reported by 
Wang et al. (Cf^-^ = 8.57Teexp[-0.82A(0)/A;BT] for T < 0.4Te). This 
discrepancy is supposed to be indicative for the opening of an additional 
gap below about 10 K . The solid line in figure 1 indicates that a fraction of 
electrons corresponding to a normal state 7 ~ 1.4 mJ/molK^ is tentatively 
accounted for by the BCS-fit with A(0)//cbTc ~ 1.9 while a second fraction 
corresponding to 7 ~ 1.0 mJ/molK^ contributes to the smaller gap opening 
at temperatures well below Tc/2. 

To further screen the quality of the sample, temperature and magnetic 
field dependent resistivity measurements p(T, H) were performed from 0.5 K 
up to room temperature. Shown in figure |^ (right axis) is p(T) at zero 
field. The transition into the SC state occurs at = 38.9 K, which is in 
fine agreement with already published data. The RRR ratio of this poly- 
crystalline sample is about 6. The resistivity behaviour in the normal state 
region matches a dependence according to p(T) = pq + AT"^ with the residual 
resistivity po = 12.5 pQcm and the coefficient A = 9 x 10~^ pQcm. A be- 
haviour of p(T) was recently reported and it seems to reflect interactions 
between charge carriers. A study of the Hall coefficient implies that electrical 
transport is dominated by holes [1^. The value of the coefficient A, however, 
is significantly smaller than that known e.g., for highly correlated electron 
systems, but seems to reflect the modest density of states at Fermi energy 



11|. It should be mentioned that other power laws with an exponent close 



to 3 were reported for sintered material of MgB2 p^ , 0. Measurements of 
the resistivity down to 0.5 K and in fields up to 12 T indicate that Hc2 is well 
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above that limit. Interestingly, the transition region becomes much broader 
with increasing fields, but different values of imprinted currents (from 5 to 
40 mA) do not change the width of the transition. 

The temperature dependent thermal conductivity A of MgB2 is shown 
in figure ^. The overall behaviour of A(T) is typical of an intermetallic 
compound where scattering on static imperfections prohibits a pronounced 
maximum to occur, as it is the case in pure and simple metals. Moreover, 
the absolute magnitude appears to be of that order, as usually found for 
intermetallics. Anomalous behaviour of A(T) in the proximity of Tc is not 
observed and a local maximum or a pronounced shoulder below Tc do not 
occur in the investigated sample. 

Generally, the total thermal conductivity of metals consists of a sum of 
an electronic contribution Ae and a lattice contribution A^: 

A = Ae + A;. (1) 

In order to separate both contributions from the total measured effect, the 
Wiedemann- Franz law is applied, assumed to be valid, at least, in simple 
metals. This model relates the electrical resistivity p with the electronic 
contribution to the thermal conductivity Ae and can be expressed as 

K(T) ^ ^, (2) 

where Lq = 2.45 x 10~* WQK~^ is the Lorenz number. 

Using equation |1| and taking into account the appropriate values of the 
electrical resistivity in the normal state region of MgB2 (compare figure |^, 
right axis) allows to split A into Ae (dashed line, figure |^) and A; (dashed- 
dotted line, figure ||). This type of analysis indicates that both contributions 
are almost equal in the entire temperature range of the normal state region 
of MgBa. 

According to Matthiessen's rule both Ag and A; are limited owing to vari- 
ous scattering processes, which can be expressed in terms of a thermal resis- 
tivity W. In the case of non-magnetic materials, the following temperature 
dependence of the electronic contribution to the total measured quantity is 



assumed to be valid [Q: 



l/Ae(T) = We{T) = We,oiT) + We,ph{T) = f^T\ (3) 
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where the subscripts (e,0) and (e,ph) refer to interactions of the conduction 
electrons with static imperfections and thermally excited phonons, respec- 
tively; a and /? are material constants. 

Equation |^ allows to detemine PVe.o and We^ph- Shown in figure ^ is the 
electronic thermal resistivity We of MgB2 displayed in the normal state region 
up to about 80 K. The solid line is a least squares fit of the data according 
to equation ^ and the dashed and the dashed-dotted hues represent Weft and 
We^ph, respectively. Thus, the deduced parameters are a = 0.55 cmK^/mW 
and f3 = 2.8 x 10~^ cm/mWK. Obviously from figure ^ the scattering of 
electrons with static imperfections of the crystal becomes dominant as the 
temperature approaches T^. 

The relative weight of Ag and A; at the SC transition temperature T^, 
as defined from the Wiedemann Franz law, also serves to determine the 
temperature dependence of Ag and subsequently of Xi below Tg. Since within 
the BCS theory. Cooper pairs do not carry heat and entropy, the scattering 
terms of equation ^ have to be modified in order to account for the decreasing 
number of unpaired electrons. 

In the SC state, the thermal resistivity can be represented as 

w: = 1/a: = wi, + wi^, = ^ + (4) 

with t = T/Tc [0. The functions f(t) and g(t) were calculated repeatedly 
and agree well with experimental findings [|TB], |T7|, |T^. In the dirty limit of 
a superconductor, the first term of equation | dominates, i.e., 1/Ag = ~ 
W^Q. In terms of the BRT theory [||] A^(t)/A^ is a universal function of t, 



dependent on the value - and the temperature dependence of the SC gap A. 

On the contrary, clean limit superconductors are dominated by the second 
term of equation |, revealing 1/Ag = ~ W^p^^. Geilikman et al. |T3 



have calcutated and tabulated g{t) which yields again a universal behaviour 
on t. Differently to scattering on imperfections, Ag(t)/Ag(T = Tg) in the 
BCS limit increases initially with decreasing values of t in spite of a rapid 
decrease of electronic excitations. A maximum occurs at t ~ 0.28 with 
Ag(t = 0.28)/Ag(T = Tg) = 2.44. Various high temperature superconductors 
are found to exhibit a maximum in A(T) below Tg and thus the origin of that 
feature is, at least partly, attributed to a significant scattering strength of 
electrons on thermally excited lattice vibrations (compare e.g. ||18|| ). 

Figure evidences that slighty above Tg, Weft exceeds We^ph by more than 
one order of magnitude and contributes at this temperature about 95 % 
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to We- This implicitely favours a description of the thermal conductivity 
of MgB2 based on scattering of electrons on impurities. Nevertheless, for 
the present analysis of the data both terms of equation ^ are considered in 
order to analyse the electronic contribution to the total thermal conductivity. 
Taking into account the functions f{t) and g{t) and the numerical values of 
a and /3 allows us to determine Ag below (compare figure ^ panel (a)). Ag 
decreases with decreasing temperature and is primarily determined by the 
BRT function f{t). The difference A — Ag represents the phonon-originated 
thermal conductivity Af in the SC state. The latter appears to be larger 
in the SC state than Ag. For temperatures T/Tg < 0.4, the lattice thermal 
conductivity of MgB2 becomes dominant. This nicely agrees with theoretical 
considerations |T^. The lattice term A^ is constrained by various scattering 



processes; among them are interactions of the phonons with electrons, point 
defects, dislocations or sheetlike faults. 

To account for observed deviations of the SC gap of MgB2 from the BCS 
theory, the experimental data from tunnel experiments on MgB2/Ag and 
MgB2/In have been used to modify the function f{t) and thus W^q. The 
smaller the value of A(T) with respect to the BCS theory, the less steep is 
the decrease of f{t) when the temperature is lowered. Since A(0) of MgB2 
as obtained from that study {2A{0) /knTc ^ 2.4, ||) is well below the BCS 
value (2A(0)//cb^c = 3.5), Ag - in this type of analysis - becomes larger in 
the SC state down to T/Tg ^ 0.2. Still, A^ > Ag. Taking Af as derived from 
the BCS-like gap, however, provides a slightly smoother crossover from the 
SC - to the normal state region of MgB2. 

Panel (b) of figure § shows the temperature dependent thermal resistivity 
= 1/Ag. Here, the data derived from the BCS model are used. Obviously, 
scattering of electrons by phonons in the SC state contributes just a fraction 
to the thermal resistivity W^, and therefore W^^ is the most significant term 
below Tg. This, of course, will not change if the actual dependence of A(T) 
is considered. 



4 Summary 

The measurement of the thermal conductivity of MgB2 reveals no pronounced 
anomaly at T = Tg and furthermore no local maximum occurs at tempera- 
tures well below Tg. Such a behaviour is most likely caused by the dominance 
of electron (hole) scattering on static imperfections present in the investi- 
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gated MgB2 sample. In fact, an analysis based on the BRT model and the 
model of Geilikman evidences that W^q is the predominant term. Beside the 
classical BCS behaviour of the SC gap, an attempt was made to incorporate 
the actual gap behaviour reported for MgB2 p. The subsequent analysis 
shows that the electronic contribution to the thermal conductivity becomes 
larger than in the former case. The application of the Wiedemann Franz law 
to the experimental data indicates that both the electronic - and the lattice 
contribution to the total thermal conductivity are almost of the same size 
over a large temperature range. 
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6 Figure captions 




Figure 1: The electronic specific heat of MgB2, Cei{T)/T versus T obtained 
by subtracting the lattice contribution, Ciat = C^^ — ■jT, where 7 ~ 2.4 
mJ/molK^. The full line indicates a BCS temperature dependence of Cei as 
explained in the text. 
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Figure 2: Temperature dependent thermal conductivity A(T) (left axis), and 
electrical resisitivity p{T) (right axis) of MgB2. The dashed and the dashed- 
dotted lines are the electronic - and the lattice contributions Ag and A^, re- 
spectively. 




Figure 3: Temperature dependent electronic thermal resistivity Wf,{T) of 
MgB2 in the normal state region. The dashed and the dashed-dotted lines 
are the contributions due to electron- imperfection - and electron- phonon scat- 
tering o s-nd We^ph, respectively. 
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Figure 4: (a): Temperature dependent thermal conductivity A(T) of MgB2 
in the SC state. The sohd and the dashcd-dottcd hncs represent electron- 
and lattice contributions XI and Af , derived from the BCS - and the modified 
model as discussed in the text, respectively, (b): Temperature dependent 
electronic thermal resistivity, W^{T) of MgB2 in the SC state derived from 
the BCS model. The dashed and the dashed-dotted fines are the contribu- 
tions due to electron-imperfection - and electron-phonon scattering W^q and 
Wlph, respectively. 
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